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Abstract 

We consider second order inflationary perturbations in the case of two scalar fields, a 
and the inflaton (p. We derive an expression for the non-Gaussianity of perturbations 
and apply the results to hybrid inflation. We isolate the contributions due to a and 
evaluate the resulting terms to show that cr-induced non-Gaussianities dominate 
over inflaton- induced non-Gaussianities when m 2 > T]H 2 , where 7] is the slow-roll 
parameter. This may provide a useful test of hybrid inflation in forthcoming CMB 
experiments. 
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1 Introduction 



With the advent of WMAP |JJ, measurements of the temperature and polarization fluctua- 
tions of the cosmic microwave background (CMB) now begin to reach the level of precision 
needed for testing various models of inflation and hence, by implication, particle physics 
at very short length scales. The spectral index of the two-point correlator of scalar pertur- 
bations is a well-known measure of the slow-roll parameters, which are directly related to 
the form of the inflaton potential. Tensor perturbations, the spectrum of which could be 
extracted from future CMB polarization and temperature fluctuation data, would yield 
independent constraints on the potential j2j |Hj - In case of multi-field inflaton models [I], 
testing is naturally more involved. For instance, the constraints imposed by the tensor 
modes are less severe j3J. Moreover, in addition to the purely adiabatic perturbations of 
the single-field inflation, there could be a small component of entropy (isocurvature) per- 
turbation (for a general treatment of the first order entropy perturbations, see 0). There 
could also be a correlation between the adiabatic and entropy perturbations which is a 
source for additional freedom in fitting the CMB data jSj. 

In addition to the spectral features of the CMB, the statistical properties of the tem- 
perature fluctuations also provide potentially important information about the origin of 
the primordial perturbations. This is an aspect that has received much less attention as 
compared to the spectral considerations - and up to now, perhaps for a good reason. 
From observations we know that the spatial distribution of the temperature fluctuations 
on the microwave sky appear to be random so that their statistics is mostly Gaussian [Zj. 
However, in many inflationary models there are several sources of non-Gaussianities; the 
self- interaction of the inflaton is one obvious example [5]. In models with more than one 
scalar field the emergence of non-Gaussianity at some level cannot usually be avoided. 
These non-Gaussianities may be small, but they could nevertheless yield important con- 
straints on models of inflation. Indeed, observations of CMB non-Gaussianity are only 
now starting to be sensitive enough to be able to distinguish between different inflation- 
ary scenarios [Oj (see also pfj] and references therein). 

Non-Gaussianity in standard single-field inflation has been considered e.g. within the 
framework of stochastic inflation ^JJ ^] . The three-point correlation function was dis- 
cussed in [13]. For Gaussian statistics, all the iV-point correlators are related to the two- 
point correlator with the odd correlators vanishing. Hence the three-point correlator of the 
scalar perturbations is the lowest order measure of non-Gaussianity of the perturbations. 
In jH] the induced cubic self-interaction terms of the inflaton were computed together with 
an estimate for the resulting non-Gaussianity, which was found to be small. Non-vacuum 
initial states and non-Gaussianity have been addressed in [TI] ITK] ITT)] . Nonstandard infla- 
tion scenarios have also been studied in the context of non-Gaussianity, such as models 
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with higher order derivative operators |17j . or models with a varying inflaton decay rate 
[TH] . However, many of these treatments (with the exception of jS]) have been lacking in 
that non-Gaussianity has been introduced by hand, often relating it to the magnitude of 
the first-order perturbation. In general a consistent treatment of the perturbation theory 
up to second order is required to find the Fourier transform of the three-point correlator, 
called the bispectrum. 

A full computation of the bispectrum of single-field inflationary models has recently 
been performed in second-order perturbation theory by Acquaviva et al. [E2]. They con- 
sidered second-order perturbations in both the metric and the energy-momentum tensor 
and, using the Einstein equation, solved the evolution of the gauge-invariant comoving 
curvature perturbation, computed up to the first order in the slow-roll parameters. This 
enabled them to find out the gauge-invariant gravitational potential bispectrum. 

Despite the occasional incompleteness of their treatment, most single-field studies 
suggest that the amount of primordial non-Gaussianity is unobservable at present (see 
however jTHj) and possibly unobservable even by the Planck satellite J9]. However, many 
models of inflation are not based on the simple single-field approach. Indeed, one of the 
most popular realization of the inflationary paradigm is hybrid inflation [2Il] which relies on 
two scalar fields, the slowly rolling inflaton ip, and an additional field a, sometimes called 
the "waterfall field" . During inflation a stays roughly constant but eventually triggers the 
end of inflation by rapidly cascading down the potential to its minimum (see e.g. [SHE])- 

Non-Gaussianity in multi-field inflation has been studied by Bernardeau and Uzan 
[221 12H] (even before them Yi and Vishniac |21j concluded that non-Gaussian effects can 
be important in multi- field models but their model lacked compelling physical motivation). 
In their case inherently non-Gaussian isocurvature perturbations are transferred to the 
adiabatic modes during inflation. The necessary ingredients for this kind of mechanism are 
a self-coupling for a transverse scalar field for generating the initial non-Gaussianity, and a 
coupling in the potential leading to a curved trajectory in the field space (and trough that 
to a mixing between isocurvature and adiabatic modes). This approach is very specific 
and is not applicable to a generic multi-field inflation. Indeed, as in the case of single-field 
inflation, the proper treatment of non-Gaussianity requires expanding metric and scalar 
field perturbations up to second order. Second order perturbation theory [23 [2EI I2Z] has 
been employed both in single-field models [HI EH] and multi- field models j2Hj, where the 
emphasis was on the connection between primordial perturbations and CMB temperature 
fluctuation at large scales. In j2H] it was applied to a study of isocurvature and adiabatic 
perturbations in multi- field models (see also [HU]), where it was found that the isocurvature 
perturbation sources the gravitational potential on long wavelengths in the second order 
even without any bend in the trajectory in the field space. We reach the same conclusion 
in our study. Second order perturbation theory has also been applied to the curvaton 
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scenario |HU 133 E3] , where the curvature perturbations are generated after the inflation; 
for a discussion on non-Gaussianity in curvaton models, see [3U |2H EH] ■ The enhancement 
of primordial perturbations after inflation has been studied in 36\, where it was found that 
even tiny non-Gaussianities produced during single-field inflation could be significantly 
enhanced by gravitational dynamics after inflation. 

This paper is organized as follows. In Section 2 we expand both the metric and energy- 
momentum tensor perturbations of two scalar fields, denoted as (p (the inflaton) and 
a, up to second order, following the single- field treatment presented in [TH]. We then 
write down the perturbed Einstein equations up to second order and construct a master 
equation for the second order metric perturbation (2) . In Section 3 we apply the master 
equation to hybrid inflation. Section 4 is devoted to constructing second order comoving 
curvature perturbation first for a general two-field case and then for hybrid inflation. The 
contribution due to the second scalar field a is isolated. In Section 5 we study the curvature 
perturbation induced by the second field in more detail. In particular the amount of non- 
Gaussianity due to a is compared to non-Gaussianity due to inflaton. Finally in Section 
6 we discuss our results and draw conclusions. 



2 Second order inflationary perturbations of two 
scalar fields 

2.1 Metric perturbations 

Let us start by assuming a spatially flat Robertson- Walker metric in conformal time 
ds 2 = g, Jil/ dx^dx v = a 2 {r){—dr 2 + dx 2 ). The components of a perturbed metric up to an 
arbitrary order can be written as 



9oo 



9oi 



9ij 



oo 1 

E^r, 

r=l 

\ r=l ' / r=l 



X 



(2) 
(3) 



where the functions (r) , \ ?/> (r) , and Xij represent the rth order perturbations of the 
metric. 

Let us perform the standard splitting into the scalar, vector, and tensor parts by 
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writing [251 EE! 



<4 



I 



(4) 
(5) 



where <9 l c<j- r) = <9 l xj° = 0, l x« = 0, xV = 0, and = didj — \5ijd k dk- When considering 



only scalar fields up to the second order we can follow Acquaviva et al. [19] and neglect 
ujf\ xf \ and xfj, thus obtaining the metric 



9oo 
9oi 

9ij 



-a(r) 2 (l + 20 (1) + (2) ), 



(6) 
(7) 



a r 



(1 _ j^W - ^) tf y + ( + ^ 



-~(^xf + ^xf + xg ) ) 



We will adopt the generalized longitudinal gauge [23 EE!, an d set 

and x (1) = X (2) = 0- This renders the metric into the form 



ujW = m (2) = 



9oo 
9oi 

9*3 



-a(r) 2 (l + 20 (1) + (2) ) 







a r 



(i - - r } ) ki + i (ftxf + ^xf + X?) 



(9) 
(10) 

(11) 



which is the basis of our subsequent analysis. The Einstein tensor resulting from the 
metric Eqs. fjH|)- ([TTJ) is presented in Appendix 1X1 



2.2 Energy- momentum perturbations of two scalar fields 

Let us now consider the energy-momentum tensor for generic two scalar fields, ip and 
a, minimally coupled to gravity and given by 37 

T^v = d^ip d u ip + d^a d v a - Q g a(3 d a p dpip + ~ g a(i d a a d p a + V(<p, a^J , (12) 

where V((p, a) is the potential for the scalar fields. In our notation, p> is the inflaton. The 
fields can be expanded up to the second order in perturbations, where the zeroth order is 
the homogeneous part denoted by the subscript 0: 



<p(r,x) = <po{t) + 6 w (p(t,x) + -6 (2) <p(t,x), 
(t(t,x) = (To(t) + 5 {1) a(T,x) + ^5 {2) cr( y T,x). 



(13) 
(14) 



Similarly, the energy-momentum tensor can be expanded up to second order as 

T^ v = T^ 0) + 5 m T" u + ^ (2) T^, (15) 

where denotes the background value. The components of the energy-momentum 

tensor (JT2*|| up to the second order are presented in Appendix iBl 

2.3 Einstein equations 

Using the results of the two previous subsections we may expand the Einstein tensor up 
to the second order as 3 



(16) 



We list the components of the Einstein tensor for the metric (JHJ) - (|lip up to the second 
order in Appendix 1X1 

The Einstein equations can now be written in the component form. The background 



equations 4 G°q 0) = k 2 T° (0) and G^p = k 2 T %i - 0) are found to be respectively given by 

.2 



;i7) 



and 



It 



H 2 + 2H' = -y (V 2 + of) + K 2 aV • 

Here Ji = a' /a and n 2 = StxGn = Mp 2 , where Gn is Newton's constant and Mp is the 
reduced Planck mass. 

The first order perturbed equations 5 W G° = k 2 5 w T° , S^G^ = k 2 5 {1) T% and S^G^ = 
K 2 5 w T l j take respectively the forms 

6H 2 <p {1) + 6Hip w ' - 2d l d i ?P (1) 

._2 



\uip oa 



2Hd i( t) w + 2<9^ (1) ' = k 2 {<p' Q «9,<5 (1 V + <?' Q d^a) 



(19) 



2H(p m ' + 4— - 2H 2 <J) W + d k d k <J) m + 4H</> (1) ' + 2^ (1) " - d k d k ^ m J 5' 



+ <9 4 <9, (^ (1) -0 (1) ) = 



2^(1) 



<9V 



<9U 



a 2 ^5 (1 V + ^S m a 



dip 



da 



3 Note that Acquaviva et al. ^Hj expand — G^l + S^G^u + <^G>„ without the factor one half 
in the second order term, although they have it in the expansion of the energy-momentum tensor. 



4 Note that G° (0) = GV = T° (0) = tT 



0. 
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where we have defined w 2 = (p' 2 + a' 2 . 

In the second order perturbed equations we use the well known |3S] first order result for 
the scalar fields and set ^ (1> = The components 5 {2) G% = k 2 5 (2) T%, 5 w & = k 2 5 (2) T 1 
and 5 {2) G l • = k 2 5 (2) T*- read respectively as 



3 3 

a" 



H 2 (f) i2) + — (2) + 3?# (2)/ - did l ^ 2) - YIU 2 (& w f -3d t <f) w d l <J) 
a 



a 

~2~ 



dV 



dV 



d 2 V 



d 2 V 



Hd 1 ^ + <9V (2) ' + ^d k d k X i(2) ' + 20 (1)/ <9> (1) + 80 (1) d> (1) ' 



K 



(20) 



-d k d k (f>^ + H^' + — (2) + ft 2 (2) - -9 fc 9V (2) + ^ (2) " 
2 a 2 

+ 2^ (2) ' + AH 2 U m ? - 8— (0 (1) ) 2 - 8H(P W (P {1)> - 3d k( f) w d k (t) m 

a 



t 1 ) A F) k A,m 



40 (1) <9 fc <9' 



J 2 2 



1(2)' 
3 



d k d k xj ) + 2<9> (1) dj<p m + A(f) w d l dj<p w 



2 (AW\ 2 



2 (<^ ( V + a^V) + 2w 



- 2 (d k 8 w (pd k 8 m (p + d k 5 w ad k 5 w a) 



d 2 V 



d 2 V 



5 w <p5 w a 



8\ 
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In writing the 00 and ij components we have made use of the equality k 2 w 2 /2 = Ti 2 — TV 
of the background equations (fP7j) and (fT5]l. together with the relation a" /a = Ti 2 + TC 
which follows directly from the definition of Ti. 



2.4 The master equation 

Here we follow the procedure outlined in [TH] for the single field case. Our purpose is 
to derive master equations for the second order metric perturbation (2) . To this end, 
we take the divergence of the iO component of the second order Einstein equation using 
the background metric Sy, that is, we operate on Eq. (|2Uj) with <9j = S^d* (recall that 
diX l{2) — 0) together with the inverse of the spatial Laplacian, denoted by A -1 . The result 
is 



where 



1 



-(cp'^cp + a'^a) 



K 



(21) 



ai)' 



did 



+ 10d^ (1) '<9> (1) + 8 6 m d i d i ^ 1) 



+ 2 ^diff (<p' Q 6 m <p + a' 5 m a) + 2 drfVff (<p' 6 w <p + a' 5 m a) . 



(22) 



Taking the trace of the ij component of the second order Einstein equation and re- 
calling that xT = d l Xi = we obtain 



Idid^ - -did 1 ^ = -n<p {2) ' - — (2) - n 2 <p (2) - ^ (2) " - m^' + 8— (0 (l) ; 

3 3a a 



- AH 2 (0 (1) ) 2 + 8H0 (1) (1) ' + Id^d'^ + l^did 1 ^ 



+ 



(0 (1) ') 2 + (vW + ^a') + \ ((« ( V) 2 + (* (1 V) 2 



2(<p' 5 m <p' + a' 5 m a')(f> m + 2w 2 



- \ (d^ipd^cp + d^ad^a) - °- 



dV „ l2 , dV 

— 5 (2 V + a 
dip oa 



i¥V (^ (1 V) 2 + S (^ (1) ^) 2 + 2 ^T<5 (1 V s m * 

do 1 i-i- 



d<p 2 



dip da 



(23) 



From Eq. ()21)1 and its derivative with respect to time we obtain an expression for the 
second order terms of the metric perturbations on the right hand side, which transforms 
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the Eq. (231 into 



-d t d l <f) {2) - ~<9;<9V (2) = 8— (0 (1) ) 2 - m 2 U w f + 8H<p m <p (1) ' + 
3 3a 3 



+ -0 (1) <9 4 <9> (1) + + A~V + 2HA- 1 a - k 2 A-^' 

3 

- 2Hk 2 /\- 1 (3 + |i ((5 ( V) 2 + (5 (1 V) 2 



- 2 (<p' 5 (1 V + 5 (1 V) <p w + 2w 2 (0 (1) ) 2 



- (di5 (1) <pd l 5 w <p + diSVa&SVo) 
6 



a 

T 



s 2 v 



(5 (1 V) 2 + (^ (1 V) 2 + 2 ^-^V^'rr 



d 2 V 



dip 2 v " r ' da 2 y " ' ~ dip do 
Taking the inverse Laplacian of the previous equation ()24j) leads to 

^ = ^ _ A -i 7 ; 



(24) 



(25) 



where 7 is three times the right hand side of Eq. (|24|). We then plug this result into Eq. 
(J21j) and obtain 



** + H^ + A->a A - V _4 A -y. (26) 



Thus far we have made use of the iO and ij components of the second order Einstein 
equations. By utilizing Eq. (|25[) together with the derivative of Eq. ()26|) with respect to 
the conformal time r the 00 component of the second order Einstein equation can be 
written as 

(2) " - did 1 ^ + 2H4>^' + 2H'4> {2) = 12H 2 (0 (1) ) 2 + 3 (0 (1) ') 2 + 80 (1) <9 4 <9> (1) 
+ 3d l <p w d i <f) m + 2HA~ 1 a - 2Hk 2 A- 1 (3 - 2HA- l i - A~ V + ^A" 1 ^' 



+ A-y + 3WA-V - 7 + k 2 5 (2 V + 0% S (2) a - - [(5 (1 V) + (S w a'Y 
- ~ (d^ipd'S^ip + d^ad^^a) - 2w 2 (0 (1) ) 2 + 2 WV + <P 



(i) 



a 

T 



( ,. v)2 + ( ,.v) 2 + 2 |^,,V,.v 



(27) 



where we have also applied the background field equations |HZ] ip'o + 2H.ip' + a 2 dV/dip = 
and (Tq + 2^ + a 2 dV/da = 0. 

Eq. (J27|) is our master equation describing the evolution of the second order metric 
perturbation (2) . We have not made any approximations in deriving it; therefore it applies 
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generally to all inflationary models with two scalar fields. Understanding the evolution of 
(2) is of key importance since it is going to be used later in calculating the second order 
gauge invariant curvature perturbation 1Z (2) . 



3 Hybrid inflation 

3.1 Basic features of the model 

Until now our treatment has been very general and we have not made any assumptions, 
e.g., about the form of the potential. In the second order we now have three different 
equations but four unknown functions (2) , ip (2) , S {2) (p, and 5 (2) a. In fact, as is noted in 
[3H], the dynamics of a system of two scalar fields can not be described by just one 
equation even in the first order. Therefore, we need an additional constraint in order to 
solve completely the dynamics of the system. 

In our case such a constraint is provided by the fact that in hybrid inflation models 
the second scalar field a, the "waterfall field", sits at the bottom of a non-curved valley so 
that we may take <To = 0. Consistent with this we also assume that the inflation potential 
does not contain any terms linear in a so that <j = is indeed a local minimum 5 . As we 
will discuss, in hybrid inflation the perturbations in a are decoupled from the evolution 
of (2) to the extent that the dynamics of 5a can be considered independently. 

The potential of the hybrid inflation scenario reads [I] 

V((p, a) = V - \my + ^Aa 4 + X -m\ 2 + \g 2 a\ 2 . (28) 

As the parameters are related by Vq = m^/AX, the number of free parameters is four. 

The critical value of the inflaton field defines the point where the "valley" in the a 
direction disappears. It is given by 



<P 1 = ^ ■ (29) 

Inflation ends when ip reaches the critical value. Depending on the potential in the a 
direction, there could be a few additional e-folds due to the slow initial motion of a. The 
usual slow-roll parameters are defined by [HZ1 Hi] 

e= J_fI^V = (30) 

2k 2 \V dtp J " H 2 2 H 2 ' K } 



k 2 V dip 2 Hip 



3 Note that our formalism would apply also to the Linde-Mukhanov model |40j . 
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where the dot denotes derivation with respect to the cosmic time, and H = a/a is the 
Hubble parameter. During inflation, i.e. when ip > (p c , o ~ and Vq dominates the hybrid 
potential. Hence, using the potential (jZH|) the slow-roll parameters read 

2 



1 

2^2 



m 2 ip 



m 



V 



(32) 



During the slow roll, the inflaton field receives perturbations with a power law spectrum 
Vtp = k 3 \5ip\ 2 /2n ~ k n ~ x . The almost scale invariant nature of the perturbations is guar- 
anteed by the smallness of the slow-roll parameters with n — 1 — 6e + 2r], where n = 1 
corresponds to the exact scale invariance. 

By virtue of the construction of the model, the perturbations generated by the inflaton 
are mainly Gaussian. Terms quadratic in ip produce perturbations proportional to v?o 5 {1) <p. 
The waterfall field a does not give rise to this kind of perturbations since its background 
value <7 = 0. Therefore the leading order perturbations are Gaussian and due to the 
inflaton. However, as we will show below, o may well induce second order non-Gaussian 
perturbations which are larger than the small second order inflaton perturbations and 
could be observationally significant. 

Focusing now on hybrid inflation and setting cr = (together with d 2 V/dadip = 0), 
Eq. (j27jl can be written as 

(2) " - di&Q™ + 2H(f) {2) ' + 2H'(P {2) = 12H 2 (0 (1) ) 2 + 3 (0 (1) ') 2 + 8</> (1) «9,# 
+ 3^0 (1) 9> (1) + 2HA- 1 a - 2Hk 2 A- 1 [3 - A~V + ^A" 1 /?' 

1 
2 



+ A-y + wA-y - 7 + «M vo' J(2 V 
i 



(5 ( V) 2 + (5 (1 V) 2 



a 2 \d 2 V ( 2 . d 2 V (1) 2 



+ 2<p' 6 w (p 



(33) 



where the definition of a does not change but (3 and 7 simplify somewhat with respect to 
their original forms (see Eqs. (J22|) . (|25|l). From Eq. (J26|) we obtain an expression for 5 {2) (f 
which we then plug into ()33|) . The result is a master equation for the second order metric 
perturbation (2) in hybrid inflation: 



1(2) « 



/-•)' + 2 ( W - 



U2) 



12H 2 , 

- A~V 



+ 3 
2k 2 



(1) ') 2 + 



Mdi&t™ + m<P w «9> (1) + 2 ( H + ^ ) a- x « 



n + ^ ] A- 1 /? + k 2 A -1 /?' - 7 + ( W - 2^ ) A-W + A" 1 / 
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1 

2 L 



(5 (1 V) 2 + (^ (1 V) 2 



> 2 uw\ 2 



2 |_fy> 



— (* (1 V) 2 + ^ (5 (1 V) 2 



(34) 



The important point is that now all the source terms are quadratic in the first order 
perturbations. Note that in Eq. (J34)) we have not made any approximations. 

3.2 Metric perturbation in hybrid inflation 

It is convenient to rewrite our master equation Eq. (|34j) in terms of cosmic time dt = adr. 
This makes it easier to make use of the slow-roll parameters (|3U|) and (|31j). We find 



i (2) + H 1 - 2 



^ i(2) • 2# 2 1 H ^0 



-2AH 1 1 + 



H<p 

H 
H~ 2 



H 2 Hip c 



!> (2) - -^<9> (2) 



) 2 - 24#0 (1) (1) - 4^0 (1) <9> (1) 



+ 2/7 3- 



2^A" 1 7 + A- 1 7 

00 



1o 

if 0c 



a a 



a a 



(5 (1 V) 2 + (8 w af 
d 2 V 



-k z |2 



+ 8^ (0 (1) ) 2 - 80 O (1) 5 (1 V 



(£ (1 V) 2 + ^ (^ (1 V) 2 



(35) 



where we have written —7 explicitly (see Eq. ()25Jl). 

We need to evaluate the term A -1 (K 2 f3/a — a/a) in Eq. (|35|). In terms of the cosmic 
time t we may simply write 



- = 2 $0* ( ) + 6 ( <9^ (1) <9> (1) + <jF>di&* 
a 



(36) 



(see Eq. ()22j) ). To derive an expression for (3 /a we make use of the first order Einstein 
equation ()19|) . Setting = we obtain from the 00 component and from the Oi 
component an equation of motion for the first order metric perturbation given by 



i« + = — 00 <y (i v 



(37) 
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After some algebra we can write 



t = Hi (5 (1 V) 2 + 3 ip di&8V>tp + 3 ip 84 m &8 w tp 
a 2 0o 



+ 



2 did k d k <j> m &6 w (p + d k d k <j) m di&SWy 



K 2 a 2 (fo 



K 2 a 2 (po 



+ diS w a &8™d + 8 (1) & d^S^a 



(3f 



Thus we obtain 



a ay 2 0o 

+ A" 1 (<9A<9 fc (1) d l 5 w if + d k d k ct) m di&SVtp) 
a 2 (p 



(39) 



where we have again utilized the equation of motion Eq. (pT7|) . 

In hybrid inflation constraint a decouples completely from the first order perturbed 
Einstein equations Eq. (|19J) and we may isolate the contributions due to a from those 
that are due to the inflaton <p. Taking into account Eqs. (|22|). ()24|1. (J25j) and (j3U|) we thus 
write 



a = a 



ft 



P = P v + P <r = p ip + ad i (SW&tfSVo-) , (40) 
7 = lv + 7<T = lv - %a 2 K 2 HA- x di (6 w a&6 m o) - 3a 2 K 2 A~ l di (S^&^a)' 



, 2 2 

+ an 



2 v 7 2a 2 2 9a 2 V 7 



and 



A- 1 ( k 2 — — — 

a a 



A- 1 Ik 2 — — — 
a a 



+ a-M^-- 

a a 



A- 1 U 



2 P _ a 
a a 



AC' 



(41) 



where the subscripts ip and a denote respectively the inflaton and a contributions. The 
inflaton contributions have all been computed previously by Acquaviva et al. ^H] , whose 
treatment is still valid here due to the decoupling of ip and a. Plugging Eq. (jUJ) into 
Eq. (J3*5j) and dropping terms next to the leading order in slow-roll parameters the master 
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equation Eq. (f33j) can be written as 





(42) 



The inflaton part can be computed from Eq. ()35|) but it has already been done in |19j . 
The last two lines in Eq. ()42j) represent the second order metric perturbation solely due 
to the second field of hybrid inflation. 

We have now succeeded in completely isolating the contribution of o to the dynamics 
of the second order metric perturbation (2) within the hybrid inflation paradigm. This is 
of utmost importance since the master equation, Eq. (}4*2*j) . plays a key role when we later 
compute the curvature perturbation. The decoupling of ip and a makes it also possible to 
isolate the contribution of a to the curvature perturbation also. 

4 Curvature perturbations 
4.1 General formula 

Let us now derive the second-order gauge-invariant curvature perturbation for two scalar 
fields following the procedure described in [THj (for an exact treatment see [12] )• We denote 
the gauge-invariant curvature perturbation by 1Z and expand it up to the second order in 
the already familiar way 6 



We are mainly interested in the second order part. Our starting point is the first order 
quantity [5] 



Instead of the first order quantities we write the expansion up to second order for the 
metric perturbation ip = tp (1) + |V (2> an d the scalar fields 5(p = 5 m <p + \5 {2) ip and 5a = 
5 w a + \5 {2) a. We obtain 




(43) 




(44) 



6 Recall that the curvature perturbation does not have a homogeneous part. 




(45) 
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Consider then the following second-order shift of the time coordinate [213 ESI ITH] 



T — > T 



(1) 2 V (1) (1) (2) / ' 
which transforms ip {2 \ 5 (2 V an d 5 {2) <j into 

^ = ^> + 2^ (</> (1) ' + 2?#«) " {W + (^j) 2 

= 5 (2 V + + ¥>o£(V + 25 ( V) + ^ 

= 5 (V + ^ 1) (a% + ^; + 2W)+^ 

Thus, the expansion ()45|) transforms as 7 



(46) 



(47) 



(2) 5 



'(2) 



/ 2 i / 2 
^0 + ^0 



R« + - 



n 



ip' 5^!f + a f Q 5^a 

/ 2 . 72 
<^ +CT 



+ ^(2) 



/ 



'2 W 1 )) 



V <^o 2 + 0", 

W + 2H 2 -H 



Wo + °o°o 
^ + < 2 



- (2fru 



where, following [T^], we have denoted 



y> a<V + ct 

^ 2 + ^ 2 



, (48) 



(49) 



By virtue of the first order transformations ip w = — 'H& 1) i <5 (1 V = 5 (1 V + ¥>o£(i) 



and tf^a = + a' f ° we find 



The expansion 



T - T = £ 
J ^ ^(i) 

can now be written as 



W + 2H 2 -H 



^ + ^ 2 



(50) 



+ a' n Sa 



*p + 
1 



V ^ 2 + ^ 2 



T + T C 



»(i) 



6 



(28^ - d%,) 



(51) 



7 The first order part IZ 1 - 1 ^ remains unchanged under the transformation Eq. 146|l . 
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We also solve £ ( 1} from Eq. (JH0|) and insert it into the T + T term above. Note that by 
virtue of the first order transformation uj {1) = uj <1) — the last term can be written as 

- ~ (2d i cu (1) - d^J = - ~ (dV 1 ' <V (1) - d&djvM) . (52) 
Therefore, after some algebra we see that 

~ / ^+^\ 1 _ 1 ~ ~ 

= ( ^ + a>6a \ 1 _ 1 } 

V fo + CT o / 

The treatment above shows that the comoving curvature perturbation 1Z = 1Z {1) + 
7j7£ (2) , which is invariant under the time shift r — > r — + ~ (£m £m — Cm) j rea ds in the 
case of two scalar fields as 



T"! — 77 



2 



where 



+ 2 H' + 2H 2 -H y °?° +g °r° 

Vo +°"o 

- -dV 1 ' <V (1) , (54) 
6 

K «.»^,» +w ^>±4^. (55) 

ip +a 

This result coincides with the one obtained in |2H] once one takes into account the field 
redefinitions there. 

4.2 Second order curvature perturbation in hybrid inflation 

With the general longitudinal gauge, here essentially = 0, and the hybrid inflation 
condition ctq = the curvature perturbation acquires the same functional form as in the 
single-field case [TH] 

where 

= ^ + n -^ . (57) 
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Rewriting 1Z (2) in terms of the cosmic time dt = adr and applying the condition = 
we obtain 



K 



H 5 ^ + r> + 



km 



+ 2H<p w + H5 m <p/& 



^ H 2 [2 + H/H - (po/Hip Q 

Making use of the relations (|25j) and (J2lj)) we find 
2H 



(5? 



K 



(2) 



LA2) 



+ H(j) (2) - A' 1 



2 (3 a 

K 

a a 



+ 



L(2) 



2H 



A _1 7- A- x 7 



+ 



km 



+ 2H<f) w + H5 w <p/<Po 



H 2 [2 + H/H-<p /H<p 



(59) 



Acquaviva et al. ^H] point out that the last term gives a subdominant contribution in the 
single field case. Moreover, it does not contain any dependence on a, not even implicitly 
through + 2H<f) {1> . Therefore, in what follows we shall neglect this term. 

Since 2H 2 / h 2 ^ — 1/e, the term (2) outside the square brackets is subdominant to the 
one inside; hence we discard it. Thus, up to the leading order in the slow-roll parameters 
we may write the curvature perturbation as 



K 



(2) 



2H 

K 2 <j) 2 



4(2) 



+ H(j) (2) - A" 1 



K 2 - - - 

a a 



2H 

K 2 <j) 2 



A^j- A~S . 



(60) 



4.3 Curvature perturbation due to a 



As we did with the master equation, Eq. (jl2j) . we may isolate the contributions coming 
from the inflaton and a also in TZ. We thus write the comoving curvature perturbation as 



K = n w + l -n {2) 



-Jim , I<£(2) , I<£(2) 



Kv + -n 



(61) 



TZ 9 contribution has already been calculated by Acquaviva et al. [15] . They plug the 
expressions (I22J) and (|25jl for a, (3, and 7 into Eq. (f6T)|) while (2) + H<p <2> is solved from 
their master equation 8 . They also take into account the fact that at large scales, k -C aH, 
can be taken constant and 



(1) 



^6^ = eH^ 
2 H * O 



(62) 



3 This is the part marked inflaton contribution in our master equation l|42|l . 
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This makes it possible to set ip {1> = eTZ {1) so that the result can be written in a deceptively 
simple looking way as [19] 



(63) 



where 



X, 



1 



1 



d^ w d^ w dt 



i) m ) dt + (e - 77) A^diK^d^ 



(64) 



The important point to stress here is that the single field contribution to the curvature 
perturbation, including the integral part Z^, is proportional to the slow-roll parameters 
and hence naturally small in hybrid inflation. 

However, in hybrid inflation the waterfall field a yields an additional contribution to 
71. We may calculate it from Eq. (|6T)|) by plugging in the a dependent parts of a, (3, 7, 
and integrating the a dependent part of our master equation (jUJ) to obtain (0 (2> + H(f) {2) ) a 
(see ^1 for the inflaton part). Leaving 7^ implicit for the moment the result reads 



n 



(2) 



2H 

K 2 ip 2 



6K 2 HA- l di(5 w ad % 5 m a) dt + An 2 A- 1 d i {5 (1) ad i 5 m a) dt 



2k 2 (5 w &) 2 dt+ [ K 2 ^(5 w af dt- [ 2^A- 1 i a dt 
J da 2 J (p 

+ J A''% dt-K 2 A- 1 d i (5 m ad i 8 {1) a) 

6k 2 HA-% (5 (1 W<5 (1) o-) +4K 2 A~ 1 d i (5 w ad l 5 w a) 

dt 



k 2 ¥o 



1 

7h 



-2k 2 (6 w &f + K 2 m 2 a (5 m af - 2^A~ l % + A" 1 ^ 



(65) 



where we have used the shorthand notation m 2 = d 2 V jda 2 and the slow-roll relation Eq. 
flrSOj) . Since 7^, Eq. is not in very transparent form, we rewrite it before substituting 
it into Eq. ([65)1 . For instance, it is not obvious whether A~ l/ j a has terms that blow up 
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outside the horizon. To study this let us write 



To 



9a 2 K 2 HA-% (5 (1 W5 (1) (t) - 3aVA-^ (5 (1) d^5 (1) a)' 



+ 3aV 



1 (§Wa) 2 - -l-diS w ad l 5 w a - (5 (1) a) 2 

2 V ' 6a 2 2 9a 2 V ; 



2 ..2 a -1 



- 3aVA 



3#d; (8 w &tf8 w o) + .9, (5 (1) o-9^ (1 V)' - ^0 (8 w a) 



I ??? ^ 

ba z 2 



2 a -1 



K A 



39i (d k d h 6 w <T&6 w <r) + -9*9* (S^V^V) 



(66) 



where we have used the equation of motion [37] 5 (1) o"+3if<5 (1) cr— a _2 (9j(9*5 (1) cr+m 2 5 (1) cr = 0. 
From the last form of Eq. we see that 7 CT indeed is suppressed outside the horizon: 
the Fourier modes are proportional to the wavenumber squared, k 2 , since the inverse of 
the spatial Laplacian, A -1 , cancels the scale dependence caused by two of the four spatial 
derivatives in each term. 

Substituting Eq. (J66|) into Eq. (J65j) we can finally spell out explicitly the contribution 
of cr to 1Z: 



K 



(2) 



Ij QHA-% (5 w ad l 5 {1) a) + AA^di (6™&&6 w o) ' 

- 2 (5 w a) 2 + ml (5 w af + (e-rfj QHA~ 2 di (<9 fc <9 fe 5 (1 W5 (1) cr) ' 
+ (e - 77) H/±T 2 di& (d k 5 w ad k S w a)' - 3A- 2 d { (d^S^ad^^a)' 

dt - A- l di (5 w ad i 5 w a) 



1 



- -A- 2 ^ (d k 5 w ad k 5 w ay 



+ 3A- 2 d, (d h d k S w <rd i 6 w <Ty + -A" 2 ^ (d k 5 m ad k 5 m a)' 
+ 3eHA- 2 di (d k d k 5 w ad l 5 w a) + ^A^dfi 1 (d k 5 w ad k 5 w a) 



(67) 



where we have used the slow-roll relations (|3U|) and (|31|1. Various complications arising 
from the derivations with respect to time and space, inverse spatial Laplacians and es- 
pecially from the integration over time still remain. The main point, however, is that we 
have got rid of terms containing the metric explicitly and every term is quadratic in the 
first order perturbation of cr. 
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5 Large cr-induced non-Gaussianities 



An exact evaluation of the cr-induced curvature perturbation Eq. (|H7jl is beyond the scope 
of the present paper. However, we can consider the orders of magnitudes and relative sizes 
of the various terms appearing in Eq. ([67)1 . Let us first focus on the first term outside 
the time integral. Since A- 1 d l (5 m &d i 5 w a) = A^d^ad^^a + A'^adi^S^a, we 
see that Eq. ()67|) actually contains a term of the form (k 2 /eH)A~ 1 di5 {1) &d l 5 (1) a. Outside 
the horizon, where k <C aH, the first order perturbation 5 W & follows the equation of 
motion 5 w a + 3H5 (1) a + m 2 5 {1> a = 0. If we make the usual assumption that the mass 
of a is small, i.e. m a < H, the non-decaying solution is proportional to exp(— m 2 t/3H). 

2 

Therefore we can estimate \5 W &\ ~ ^f-\5 (1) a\ so that 



eH 



-A^d^ad^^a 



H 2 



,6 m a 

00 



d l H 



6 m 



a 



(6* 



where we have used the relation e = K 2 ifl/2H 2 . The first order perturbation of an effec- 
tively massless field a is the same as the first order perturbation of the inflaton field <p 
during horizon exit. Hence outside the horizon we may estimate 



H Sa ° 










00 




00 



\1l 



cm 



(69) 



As a consequence, we finally obtain an estimate for the simplest non-integral term in Eq. 
(EZD: 



K 

7h' 



■A- 1 d i S w ad i 6 w a 



H 2 



A djK (1) d l K 



(!) 



(70) 



This result can directly be compared with the last term in in Eq. (|64p. Because the 
spatial derivative operators 9 in A~ x d(R, w d l lZ {1) cancel each others scale dependence, we 
can approximate roughly | A~ 1 diTZ w d l TZ w \ ~ |7£ (1) | 2 . This coincides with the statement in 
Acquaviva et al. [T^j that X v is of the same order in slow-roll parameters as the contribution 
of the (K m f term to Tlf . 

In a similar fashion we can approximate all the other terms in TZ^ which are not inte- 
grated over time. Inspection shows that in every non-integrated term the spatial derivative 
operators cancel each others contribution to the scale dependence. Replacing also time 

2 

derivatives with the factor ^ we see that there are only two kinds of terms, those pro- 
portional to |7£ (1) | 2 and those proportional to e |7£ (1) | 2 . 

The estimation of the integral part in Eq. (jfi7j) is more involved. Because of the inte- 
gration over time the Hubble parameter H, slow-roll parameters e and i], mass m a and 
perturbation 5 w a can no longer be assumed to be constants. Nevertheless, we may argue 

9 That is, the two derivatives di and the inverse of the spatial Laplacian A -1 . 
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that the evolution of these parameters is relatively slow. In order to get an estimate for 
the integral part without too many complications we take the time integration to start 
after the scales in question have exited the horizon. We also assume that the quantities 
H, e, 77, m a and 5 w a change slowly and, for our estimation purposes, take them to be 
constants. The error made this way is not likely to be significant. 

With these assumptions the terms in the integrand with an overall time derivative (or 
two time derivatives) are tractable. They produce contributions proportional to jjfc \TZ m \ 2 
and (e - rj) \TZ W \. We are then left with three terms in the integrand, which we rewrite 
using the same technique as with 7^ in Eq. We obtain 

QHA^di (5 w ad l 5 w a) - 2 (5 (1) af + m 2 a (5 w a) 2 
= A' 1 \6Hdi (SW&tfSWrf+mlditf (S w a) 2 } - 2 (5 w af 
= -2A~% (5 w ad l S w a) - 2 (S w a) 2 , (71) 

where we have used the equation of motion 10 5 m a + 3H5 (1) a + m 2 a 5 m a = 0. These two 
remaining terms in the integrand we estimate like the terms outside the integral. They 
both have two time derivatives and no overall scale dependence, so we can estimate them 
roughly as {m 2 a /H) 2 |<5 (1) <r| 2 . Hence we may write 



eH 



J ^HA^d, (<5 (1 W5 (1) a) - 2 (5 W &) 2 + ml (5 (1) a) 



dt 



4 

^ AN ^ \TZ (1) \ 2 , (72) 

where AiV denotes the number of e-folds that the scale in question has spent outside 
the horizon. The relation AiV = J Hdt follows directly from the definition iV = ln(ae/a) 
where the subscript e denotes the end of inflation. Since the maximum value of AiV ~ 60 
|3*7] we can rather safely say that the contribution from Eq. (JT2*|) does not dominate over 

m 2 a iT->m|2 

H 2 



Hence, combining the contributions of both fields, ip and a, we may write the final 
result for the order of magnitude of non-Gaussianity in hybrid inflation as 



2 

m 



n w = [ arj + be + Cal JL)(Kwy , (73) 



10 We have omitted the spatial derivative term a~ 2 did 1 a from the equation of motion. It may give 
some contribution immediately after the horizon exit but its total effect in the integration over time is 
negligible. 
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where a, b and c CT are constants. Both fields (p and a contribute to constants a and b but 
c a is solely due to the a-field. In principle, the constants could be evaluated using Eqs. 

(jEH) ) and (jHZI). Simple estimation suggests they are of the same order. We should 
however emphasize that Eq. (fT3j) should not be taken to suggest that non-Gaussianity 
in hybrid inflation is x 2 -distributed; in fact, the leading terms are non-local, as discussed 
above. 

Thus we may conclude that, barring accidental cancellations in 1Z {2 \ the second scalar 
field a of hybrid inflation gives the dominant contribution to 7Z {2> if the mass of the a 
field is large enough compared to the slow-roll parameters, i.e. if 

2 

711 

^>0(e, V ). (74) 
Since usually it is assumed that e 77 the condition (f74"j) can be written as 

ml > V H 2 . (75) 

The approach followed in this Section does not apply if a is too massive. When m a > 
3H/2, the a perturbations acquire a maximum scale dependence fc 3 |5 (1) cr| 2 oc (k/aH) 3 . In 
addition, the amplitude of perturbations |<5 (1) cr| 2 becomes suppressed by m a /H jJT]. The 
strong scale dependence is due to the fast decay of the perturbation outside the horizon, 
which can easily be seen from the equation of motion whose solutions decay proportionally 
to exp(— 3Ht/2) when m a > 3H/2. Then immediately after horizon exit a perturbations 
would be suppressed as |<5 (1) a| 2 ~ (H/m a ) \5 {1) (p\ 2 . Thereafter they would decay much 
faster than the perturbations of the effectively massless inflaton field. This would most 
likely reduce their contribution to 1Z (2) insignificant. 

6 Discussion 

We have presented an analysis of non-Gaussianity in hybrid inflation using a second order 
perturbation expansion of the metric and energy-momentum tensor 11 . In hybrid inflation 
both the inflaton and the "waterfall field" a give rise to non-Gaussian fluctuations. We 
showed that the two sources can be separated and computed the comoving curvature 
perturbation 1Z up to second order, paying particular attention to the contribution of 
a to 1Z. We found that it does not affect the first order perturbation but is present 
at second order. Our main observation is that while inflaton-induced non-Gaussianities 

11 Our formalism is applicable to those multi-field inflation models for which the trajectory in the field 
space is not curved, so that one can set the background value ao = (the absence of a linear term in a 
is also required); in addition to hybrid inflation, the Linde-Mukhanov model |4L)| is another example. 
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always are proportional to the slow-roll parameters, and therefore usually small, non- 
Gaussianities induced by o do not have the slow-roll parameter dependence but rather 
scale like m 2 a /H 2 . Hence the waterfall field a is the main source of non-Gaussianity in 
hybrid inflation whenever m 2 > rjH 2 . 

In hybrid inflation H is highly adjustable; difficulties start emerging only when the 
energy scale during inflation is lowered to ~TeV [43J. Therefore no clear-cut conclusion 
can be drawn for a generic model of hybrid inflation. Note however that WMAP data 
for the spectral index n = 0.99 ± 0.04 [H] implies that 77 ~ 0.01 through n — 1 ~ 2t] 
(assuming e 77). For illustrative purposes we could adopt 77 H 2 ~ (10 9 GeV) 2 as a typical 
value for hybrid inflation. This implies that a-induced non-Gaussianities can dominate 
already for relatively low m a . Moreover, potentially they can be large. Therefore, the 
future detection of CMB non-Gaussianities could provide a powerful tool for testing hybrid 
inflation models. 

Comparison with data is however not straightforward. This is true for any proper 
treatment of the second-order perturbations, not just for the case of hybrid inflation. A 
fundamental problem is that in the presence of non-Gaussianities the iV-point correlators 
need not be related - in principle non-Gaussianity implies an infinite number of degrees of 
freedom. Therefore a generic parametrization and observational testing of non-Gaussianity 
is difficult (for a discussion, see e.g. [HI]). 

A measure of non-Gaussianity often used in the literature is the non-linearity param- 
eter / NL defined by [TU] $(x) = & G (x) + / NL [$q(x) - (<3> G (:r))], where $ is curvature 
perturbation proportional to TZ and $g oc 7t m is the Gaussian part; the angle brackets 
denote statistical ensemble average. By construction this parametrization is ideal for x 2 
non-Gaussianity. Although the a contribution to 7Z, Eq. (|67jl. is of the form first order 
perturbation squared, it does not contain an explicit factor (TZ (1) ) 2 and is therefore not x 2 
distributed. Only the first term in the inflaton contribution, Eq. (|63|h is directly amenable 
to the parametrization above (and would imply /nl oc 77 — 3e). However, the rest of the in- 
flaton contribution, Eq. (jHlJ), together with all of the a contribution, Eq. (JoTj) . contribute 
non-locally as they depend on terms integrated over time or space. 

A more appropriate parameter would be the non-Gaussianity kernel /C, which is di- 
rectly related to the curvature perturbation bispectrum as we will show below. K, is also 
used by Acquaviva et al. (THj. It could be called a "momentum dependent" non-linearity 
parameter and is defined in momentum space by [2H] (see also 19J) 



where the constant is such that ($(&)) = 0. With this definition the bispectrum for the 




d'h d 3 k 2 5 C3) (fe! + k 2 - k)K,{k u fc 2 )$ G (fc 1 )$ G (fc 2 ) 



+ constant , 



(76) 
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primordial curvature perturbation $ reads |2*8*| ITH] 



($(fca)$(fe 2 )$(fe 3 )) = (2vr) 



3 5 (3) (fe! + fc 2 + fc 3 ) [2 /C(fei, fe 2 ) P$(fc 2 ) 
+ cyclic] , 



(77) 



where Vq>{k) is the power spectrum for the primordial curvature perturbation. The differ- 
ent contributions in Eq. (|77)l depend on the overall factors in Eq. (|73)l : hence the waterfall 
field a would give the leading contribution to the non-Gaussianity kernel whenever the 
condition Eq. (J75Jl is met. 

Although the non-Gaussianity kernel JC is not identical to the non-linearity parameter 
/nl, w e may nevertheless deduce some rough observational constraints on / NL that can 
be used to evaluate the significance of the non-Gaussianities discussed in this paper. For 
instance, since WMAP [7 a established the bounds —58 < /nl < 134 with 95% confidence, 
non-Gaussianity given by Eq. (J73|) is beyond the present observational accuracy provided 
m a < H . Given the cosmic variance, detector noise, and foreground sources, WMAP is 
eventually expected to detect |/nl| with an accuracy of 20, while Planck has the projected 
accuracy of 5 9J. An ideal experiment has a sensitivity of 3 for |/nl| Hence there is a 
possibility that Planck will detect non-Gaussianities if m a /H is rather close to unity and 
the constant c a in Eq. (|73|) is larger that unity. Thus a more accurate determination of 
the constants a, b and c a in Eq. (|73j) is clearly desirable. 

We should point out that post-inflationary effects provide another source of uncertainty 
for a comparison between data and theoretical estimates. Our formalism, and especially 
the results, Eqs. (JHHj), (joTj) and (|73j). apply only during inflation. Likewise, the definitions 
for JC and /nl discussed here are for the primordial values, i.e., for the values before the 
end of inflation and reheating. We have not considered the effects of reheating and the 
later evolution of these perturbations. It is known that in multi-field models the possibility 
for generating large non-Gaussianities after inflation exists; this can happen e.g. in the 
curvaton scenario (see e.g. jSH E2] for the model and I2HJ G23 for a discussion on 
non-Gaussianity). In such a case one should take the perturbations derived here as initial 
conditions for the subsequent evolution, and follow the treatment presented in [36J and 
applied in jEHJ. 
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Appendix 



A Components of the Einstein tensor 

We take the components of the Einstein tensor, Eq. (jlfij) . from Acquaviva et al. [19J, 
where they are presented in a general gauge. For convenience we cite them here in the 
generalized longitudinal gauge for the metric (|9J) - (jllj) . 
The background components are 

(78) 



r<o(o) 

Lr 


= -A 
a 2 


(9 


G i(0) 


1 


a 


q0(0) 


_ r yi< -°'> - 


= 0. 



The first order components are 



5 (1) G° 



j« + 6-^ (1) ' - 2^9V (1) 
a 



-2-d4 m - 2<9^ (1) ' 



a 
a' 



> (1) + + 4-^ (1) ' 

a 



+ 2^ (1) " - d k d k ip m J - PdrfM + 
The second order components are 

3 (Sl) 0« + 3 ^«' _ d . d i^) _ 12 
\ a / a 

42-0 (1 V (1) ' - 3<9^ (1) 9V (1) - 8V> (1) 0idty (1) + 12-^ (1 V (1) ' - 3 (^ (1) ') 2 



5 (2) G° = 2a- 2 



/\ 2 



a 



5 {2) G\ = 2a- 2 ( -<9> (2) + (9V (2) ' + -d k d k X i{2) ' - 4-0 (1) <9 i (1) 
\ a 4 a 

+ 4^^(D^0(i) _ 2^(1)^0(1) + 4^ (1)/ (9V (1) + 8^ (1) (9V (1) ' 
a 

(^G , = 2 a~ 2 f--^0 (2) - d^' - -d k d k X f + 8-0 (1) ^0 (1) 
\ a 4 a 



(79) 
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+ 40 (1) <9^ (1)/ + 2tfj w 'd i( f) w - 4</> (1) 'c^ (1) - 4</> (1) <9^ (1) ') 
5 (2) G\ = 2a~ 2 



\d k d k ^ 2) + -0 (2) ' + 2— (2) - ( - ] (2) - -d k d k ^ 2) 



a 



a V a 



a 



a 



a 

.2 „a' 
a 

+ 2^ (1) <9 fc <9 fe (1) 



- <9 fc (1) <9 fe (1) - 20 (1) <9 fc <9V 

- 2d k ^ (1) d k tP w - ^ w d k d k tP (1> + (> (1) ') Z + 8-^ (1 V (1) ' + 4^ (1 V (1) " 



^ - ^-0 (2) + \d%^ - -d k d k X im 



J 



+ tf^drf™ + 24> il) d i d j 4>^ - 2if) (1) d i d j 4> (1) - av (1) o>^ (1) 

- d^^djipv + 39V (1) o' i V' (1) + 4?/> (1) <9 i o^ (1) 



B Energy- momentum perturbations of two scalar 
fields 

We consider here the energy-momentum tensor T^ v for two scalar fields (p and a minimally 
coupled to gravity which can be obtained in a straightforward way by using Eqs. (JEJ), 
l|T5|) . (dl|) and ©-(HD). The results are (V Q = V{<p ,a )) 
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( 1 /2 1/2 2t/ > 
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d 2 V 



i 

5 W2* 



o , 

a- 2 (-<p' d i 6 w <p-a' d i 6 w a) , 



(82) 



5 (2) T° = 2 a- 2 

+ 2(<p' d i 6 w <p + o' d i 6 w a) 
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1 

5 m T\ 
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